Doping dependence of spin dynamics of drifting electrons in GaAs bulks 
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We study the effect of the impurity density on lifetimes and relaxation lengths of electron spins in 
the presence of a static electric field in a n-type GaAs bulk. The transport of electrons and the spin 
dynamics are simulated by using a semiclassical Monte Carlo approach, which takes into account the 
intravalley scattering mechanisms of warm electrons in the semiconductor material. Spin relaxation 
is considered through the D'yakonov-Perel mechanism, which is the dominant mechanism in III- 
V semiconductors. The evolution of spin polarization is analyzed by computing the lifetimes and 
depolarization lengths as a function of the doping density in the range 10 13 -f- 5 ■ 10 16 cm , for 
different values of the amplitude of the static electric field (0.1 -f- 1.0 kV/cm). We find an increase of 
the electron spin lifetime as a function of the doping density, more evident for lattice temperatures 
lower than 150 K. Moreover, at very low intensities of the driving field, the spin depolarization length 
shows a nonmonotonic behaviour with the density. At the room temperature, the spin lifetimes and 
depolarization lengths are nearly independent on the doping density. The underlying physics is 
analyzed. 

PACS numbers: 71.70.Ej,72.25.Dc,72.25.Rb 



I. INTRODUCTION 

Spin dynamics is one of the central focuses of semicon- 
ductor spintronics. In fact, in order to make spintronics 
an usable technology, spin signal is required to sustain 
enough long time and/or enough long distance during 
transport. This is necessary to control and detect the 
spin polarization in logic operations, communication and 
storage of information. The loss of spin polarization be- 
fore, during and after the necessary operations is a cru- 
cial point in spin device design. Therefore, a complete 
understanding of the effect of temperature, driving field 
amplitude and doping density on the spin dephasing is 
essential [l|, |2j. 

Until now, the experimental investigation of the doping 
density influence on the ultrafast spin dynamics in bulk 
semiconductors have been performed at low temperature 
(T < 80 K) 0-0]. In the detailed work of Dzhioev et 
al. the dependence of the spin lifetime on the donor 
concentration at very low temperatures (T < 5 K) shows 
a very unusual behavior, characterized by the presence 
of two maxima, ascribed to the predominance of one of 
the three different spin-relaxation mechanisms: hyper- 
fine interaction, anisotropic exchange interaction, and 
D'yakonov-Perel (DP) mechanism. Recently, Romer et 
al. @ have measured the electron-spin relaxation in bulk 
GaAs for doping densities close to the metal-to-insulator 
transition, finding that at temperatures higher than 30 
K and densities lower than 8.8 ■ 10 16 cm -3 , all measure- 
ments are consistent with DP spin relaxation of free elec- 
trons since all electrons are delocalized and other spin- 
relaxation mechanisms can be neglected. From a theo- 



retical point of view, by using a fully microscopic kinetic 
spin Bloch equation approach, Jiang and Wu have pre- 
dicted a nonmonotonic dependence of the spin relaxation 
time on the donor concentration, showing that the maxi- 
mum spin relaxation time occurs at the crossing between 
the degenerate regime and the non-degenerate one [8J, |9| . 

In this work, we study the impurity density effect on 
the fast process of relaxation of non equilibrium elec- 
tron spin polarization in GaAs bulks, by using a semi- 
classical Monte Carlo technique to solve the Boltzmann 
equation 10]. We analyze the spin depolarization of 
drifting electrons at different lattice temperatures T by 
considering only the DP mechanism, which is dominant 
in n-type III- V semiconductors [J, [Tl| . This mechanism, 
effective in the intervals between the collisions, is related 
to the spin-orbit splitting of the conduction band in non- 
centrosymmetric semiconductors like GaAs [ljj . 



II. SPIN DYNAMICS MODEL AND MONTE 
CARLO APPROACH 

In a semiclassical formalism, the term of the single- 
electron Hamiltonian which accounts for the spin-orbit 
interaction can be written as 



H. 



so 



-a -CI. 



(1) 



It represents the energy of electron spins precessing 
around an effective magnetic field [B — hfl/fiBg] with 
frequency Q, which depends on the orientation of the 
electron momentum vector with respect to the crystal 
The quantum-mechanical description of elec- 
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tron spin evolution is equivalent to that of a classical 
momentum S experiencing the effective magnetic field, 



2 



as described by the equation of motion 



dt 



(2) 



Every scattering event reorients the direction of the 
precession axis, making the orientation of the effective 
magnetic field B (that strongly depends on k) random 
and trajectory-dependent, thus leading to spin dephas- 
ingdStni- 

The Monte Carlo code used here follows the procedure 
described in Ref. The spin polarization vector has 

been incorporated into the algorithm as an additional pa- 
rameter and calculated for each free carrier, by following 
the procedure described in Ref. LOj with the difference 
that in the present paper we assume the spin-orbit cou- 
pling coefficient in T- valley (/3r) equal to 8.2 eV ■ A , as 
used in Ref. to obtain a better fit with the experimen- 
tal work of Kikkawa and Awschalom [3| ■ 

In our simulations we use a temporal step At of 10 
fs and a 5 • 10 4 electron ensemble to collect spin statis- 
tics. To achieve the steady-state transport regime, we 
run the simulation code for a transient time (typically 
10 4 time steps). After that, all the spin electrons are ini- 
tialized, the spin relaxation begins and we collect data. 
All simulations are performed in a n-type GaAs bulk 
with a free electrons concentration varying into the range 
10 13 -7- 5 • 10 16 cm -3 (non-degenerate regime) by assum- 
ing that all donors are ionized. By an exponential fit- 
ting of the decay of the spin polarization we estimate the 
spin lifetime r and the spin depolarization length L [To| . 
These parameters satisfy the relation L = Vd ■ t, where 
Vd is the average drift velocity. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

In Fig. [TJ we show the electron spin lifetime r as a 
function of the doping density N, for different values of 
applied electric field F, namely 0.1, 0.5 and 1.0 kV/cm. 
In each panel, we show four curves at the following values 
of lattice temperature T; 40, 77, 150 and 300 K. 

Up to T = 150 K, the electron spin lifetime is nearly 
independent on N until N ~ 10 15 cm~ 3 , then r increases 
with the doping density. For N > 10 15 cmT z and for each 
value of the applied field, the longest value of r is ob- 
tained at T = 77 K (t M ax ~ 14 ns at F = 0.1 kV/cm). 
At the room temperature (300 K), the spin lifetimes are 
almost insensitive to the impurity density. In the investi- 
gated range of N, the system is non-degenerate, i.e. the 
electron plasma temperature is much greater than the 
Fermi's temperature (T e ^> Tp). Hence the inhomoge- 
neous broadening <| Q(fc) | 2 > is little sensitive to N, 
while the momentum scattering rate t^ 1 is proportional 
to a linear function of N @|. So, in accordance with the 
DP classical relation r oc (<| f2(fc) | 2 > r^) -1 , for high 
values of N the spin lifetime r increases with the doping 
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FIG. 1: Electron spin lifetimes r as a function of the doping 
density at different amplitudes of the applied electric field F 
0.1 kV/cm (a), 0.5 kV/cm (b) and 1.0 kV/cm (c), and four 
different values of lattice temperature, namely T = 40, 77, 
150, 300 K. 



density [llj. Moreover, for all the investigated intensities 
of the driving field, the relaxation time r has a nonmono- 
tonic behavior as a function of the temperature. 

In Fig. [51 we show the electron spin depolarization 
length L as a function of the doping density N at the 
same values of applied electric fields and lattice temper- 
atures used in Fig. [T] In particular, in panel (a), i.e. 
for F = 0.1 kV/cm, for T < 150 K, L appear to be 
a nonmonotonic function of TV, by showing a minimum 
at N ~ 2 • 10 14 cm -3 . At higher temperatures, L is 
nearly independent on the doping density. For higher 
amplitudes of the electric field (panels (b) and (c)), up 
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to T = 77 K, L is insensitive to both the temperature 
and the doping density until N ~ 10 15 cm" 3 and slightly 
increasing for higher values of N. For T 150 K the 
effect of the doping density is marginal. To understand 
the behaviour of L as a function of N it is necessary to 
consider the interplay between r and v,i in the relation 
L — Vd-T. In fact, in the investigated range of TV, the spin 
lifetime r always increases with N; on the contrary Vd is 
a decreasing function of N. The nonmonotonic behav- 
ior of L, observed at F = 0.1 kV/cm and T in the range 
40-h77 K, arises from the fact that for 10 13 < N < 2 • 10 14 
c?ti" 3 , Vd decreases more rapidly than r increases. Vicev- 
ersa, for N > 2 • 10 14 cm" 3 , t increases more quickly and 
hence L increases too. 



In conclusion, we studied the effect of the doping den- 
sity on the ultrafast spin dynamics during drift transport 
in a GaAs bulk below the metal-to-insulator transition. 
For lattice temperatures T ^ 150 K, the electron spin 
lifetime is an increasing function of the doping density, 
in accordance with the results reported in other works 
using different theoretical approaches @, Q and with the 
recent experimental results of Romer et al. @. Moreover, 
for very low intensities of the driving field, the spin depo- 
larization length shows a nonmonotonic behaviour with 
the density. At the room temperature, the spin relax- 
ation tends to be insensitive to the donor concentration. 



10" 
N [cm" 3 ] 



IV. ACKNOWLEDGMENTS 



FIG. 2: Electron spin depolarization length L as a func- 
tion of the doping density at different amplitudes of the ap- 
plied electric field F (a) 0.1 kV/cm, (b) 0.5 kV/cm and (c) 
1.0 kV/cm, and four different values of lattice temperature, 
namely T = 40, 77, 150, 300 K. 



This work was partially supported by MIUR and 
CNISM-INFM. The authors acknowledge CASPUR for 
support by the standard HPC grant 2010. 



[1] M.W. Wu, J.H. Jiang and M.Q. Weng, Physics Reports 

493, 61 (2010). 
[2] M.E. Flatte, Nature 462, 419 (2009). 
[3] J.M. Kikkawa and D.D. Awschalom, Phys. Rev. Lett. 80, 

4313 (1998). 

[4] M. Furis, D.L. Smith, S.A. Crooker and J.L. Reno, Appl. 

Phys. Lett. 89, 102102 (2006). 
[5] R.I. Dzhioev, K.V. Kavokin, V.L. Korenev, M.V. Lazarev 

et al., Phys. Rev. B 66, 245204 (2002). 
[6] M. Romer, H. Bernien, G. Miiller, D. Schuh, J. Hiibner 

and M. Oestreich, Phys. Rev. B 81, 075216 (2010). 
[7] M. Kraufi, H.C. Schneider, R. Bratschitsch, Z. Chen and 



S.T. Cundiff, Phys. Rev. B 81, 035213 (2010). 
[8] J.H. Jiang and M.W. Wu, Phys. Rev. B 79, 125206 
(2009). 

[9] K. Shen, Chin. Phys. Lett. 26, 067201 (2009). 
[10] S. Spezia, D. Persano Adorno, N. Pizzolato and B. Spag- 

nolo, Acta Physica Polomca B 41, 1171 (2010). 
[11] N.I. Dyakonov and V.I. Perel, Sov. Phys. Solid State 13, 

3023 (1971); M.I. D'yakonov, Physica E 35, 246 (2006). 
[12] D. Persano Adorno, M. Zarcone and G. Ferrante, Laser 

Phys. 66, 310 (2000). 



